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Abstract 

We perform self-consistent Schrodinger-Poisson calculations with exchange and correlation cor- 
rections to determine the electron/hole gas in a radial hetero-junction formed in a modulation 
doped GaAs/AlGaAs core- multi-shell nanowire (CSNW) which is n-/p-doped. Realistic composi- 
tion and geometry are mapped on an symmetry compliant two-dimensional grid, and the inver- 
sion/accumulation layers are obtained assuming mid-gap Fermi energy pinning at the surface. We 
show that the electron and hole gases can be tuned to different localizations and symmetries inside 
the core as a function of the doping level. Contrary to planar hetero-j unctions, conduction elec- 
, trons do not form a uniform 2D electron gas (2DEG) localized at the GaAs/AlGaAs interface, but 

rather show a transition between i) an isotropic, cylindrical distribution deep in the GaAs core (low 
0^ ■ doping), and ii) a set of six tunnel-coupled quasi- ID channels at the edges of the interface (high 

o ; 

doping). Holes, on the other hand, are much more localized at the GaAs/AlGaAs interface and 
form either i) six separated planar 2DEGs at the GaAs/AlGaAs interfaces (low doping), ii) a quasi 

•rH , 

uniform six-fold bent 2DEG (intermediate doping), or iii) six tunnel-coupled quasi-lD channels 
at the edges (high doping). We also simulate the electron/hole gas in a CSNW-based field effect 
transistor. The field generated by a back-gate may easily deform the electron or hole gas, breaking 
the six-fold symmetry. Single 2DEGs at one interface or multiple quasi-lD channels are shown to 
form as a function of voltage intensity, polarity, and carrier type. 

PACS numbers: 73.21.Fg, 73.21.Hb, 03.65.Ge 
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FIG. 1: Schematics of the CSNW geometry. The free-carrier gas is formed in the GaAs core, 
eventually at the heterointerface with the AlGaAs shell. The shell is either n- or p-doped in the 
middle with a 10 nm wide layer, and covered by a GaAs capping layer. The two thin lines crossing 
the hexagonal section represent the facet-to-facet and edge-to-edge directions on which the charge 
density is integrated in order to obtain the equivalent sheet charge density [see text and Eq. 
Starting from the center and following the facet-to-facet direction, the core has a radius of 40 nm, 
while the thicknesses of the internal and capping layers are 50 nm and 10 nm, respectively. 

I. INTRODUCTION 

Semiconductor nanowires (NWs) are rapidly evolving into functional nano-materials^. 
One key direction is the demonstration of core- multi-shell NWs (CSNWs), multi-layered 
materials where free-standing NWs, grown vertically on top of a semiconductor surface, 
are used as a substrate for the radial overgrowth of multilayers,- - -, strongly motivated by 
applications in energy harvesting^ - - and electro-optical devices^. The resulting system is 
a radial hetero-structure, similar to 2D systems, like an hetero-j unction (as in Fig. [1]) or 
a quantum well, but grown radially, and wrapped around the core. Bridging between the 
self-assembling, bottom-up, quasi-lD nature of NWs and the flexible engineering of planar 
2D hetero-structures, which opened up decades of astonishing developments in fundamental 
physics and innovative applications, CSNWs have the potential to bring the field of NW- 
based devices to the level of a new key nano-technology. 

Critical steps have been taken in this direction. Long single- crystal, defect-free, untapered 
GaAs NWs, the material class of choice for quantum transport and high mobility, have been 
recently grownii^. Epitaxial GaAs/AlGaAs interfaces and multi-layered structures can be 
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realized with high-quality interfaces^ and selective radial doping has been demonstrated^. 
Self-catalytic or catalyst-free growth protocols have been developed, thus avoiding use of 
Au, a source of deep trap formation- High- mobility and ballistic transport, similarly to 
GaAs planar structures, are therefore realistic targets in CSNWs in the short run. 

A typical geometry, as the one investigated in this paper, is shown in Fig. [TJ GaAs NWs 
grow as hexagonal crystals along the [111] direction, exposing the six {110} facets. The 
GaAs core is overgrown by an epitaxial AlGaAs shell, possibly including a doping layer, 
and by a GaAs capping layer. Surface states easily deplete the outer layers, and a 2DEG 
may localize at the inner GaAs/AlGaAs hetero-interface, which consists of six planar facets 
several tens of nm wide, interrupted by six edges. The confined electronic system has the 
six-fold symmetry of the NW used substrate. 

It should be noted that in a typical hetero-j unction the electron gas extends into the 
GaAs layer by tens of nm, due to the balance between kinetic and Coulomb energy scales. 
This is comparable to typical NWs diameters, and the electron gas on opposite sides of the 
GaAs core may strongly couple.- 4 Therefore, the free electron gas formation in a CSNW is a 
genuinely two-dimensional problem (we assume the CSNW to be translationally invariant in 
one dimension). For example, it is easy to infer that for large doping/size of the structure, 
Coulomb energy will induce charge to localize along six quasi-lD channels at the edges of 
the interface, to maximize the average electron-electron distance. 

Ab-initio atomistic methods are computationally intensive and therefore limited to very 
thin NWs, 1-2 nm in diameter. Nonetheless, they can provide important informations on 
the band-alignment and dopant properties.— _ — Modelling of larger CSNWs is also limited. 
Recently, the single-particle properties of a cylindrical electron gas have been simulated in- 
cluding strong homogeneous and inhomogeneous magnetic fields^ and k ■ p effects. 1 ^ Due 
to the prismatic geometry of CSNWs, however, cylindrical symmetry is broken, and local- 
ization at the edges of the hetero-interface may take place.— ^ Therefore, in NWs with 
diameters in the tens of nm range, the electronic system deviates substantially from the 
idealized cylindrical shape, and the electron gas may behave as a set of quasi-lD systems, 
rather than a uniform electron gas wrapped around the core. The above simulations neglect 
electron-electron interactions, though, and the confinement at the GaAs/AlGaAs interface 
is assumed, as, e.g., in a narrow radial quantum well. As discussed above, however, in 
doped structure the in-plane directions cannot be disentangled and full calculations need to 
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be performed where the confinement is self-consistently included. 22 

In this paper we discuss the formation of the electron and hole gas in GaAs/AlGaAs core- 
shell NWs which are n- or p-doped in the AlGaAs barriers. We will show that remote charge 
density which forms at the inner GaAs/AlGaAs interface may be tuned with doping and/or 
external gates changing the symmetry and localization of the distribution. The electron/hole 
gas has a mixed dimensionality, which is locally ID or 2D. Remarkable difference, from this 
point of view, is found between electron and hole gases. We finally show that the electron 
or hole gas can be easily reshaped and tuned between 2D and ID channels, which can be 
individually depleted, by a transverse electric field applied by a back-gate. 

The following of this paper is organized as follows. In the next Sec. |TTJ the physical and 
numerical modelling is sketched, and the self-consistent procedure is described. In Sec. IIII[ 
the results of our calculations are illustrated. First, we introduce the prototypical structure 
that we investigate and the relevant simulative parameter. In Subsections IIIIBI and IIII CI 
we report our results for n-doped and p-doped systems, respectively. The origin of the 
localization patterns will be discussed in Subsection IIII Dl while in Subsection IIII El will deal 
with a CSNW subject to an external electric field in a field effect transistor configuration. 
Finally, in Sec. IIVI we draw our conclusions. 

II. METHOD 

The free electron/hole gas of the CSNW is obtained within an envelope-function approach 
in a single-band approximation, including carrier-carrier interaction at a mean-field level by 
solving the Poisson equation for the free charge distribution plus the charge due to the 
remote doping. To solve the coupled system of Schrodinger and Poisson equations we adopt 
the usual self-consistent approach. Assuming translational invariance along the NW growth 
axis z, we need to solve numerically a 2D problem in the x, y plane. At each iteration, we 
first solve the effective-mass Schrodinger equation on a 2D hexagonal domain, representing 
the cross-section of the CSNW, for both electron and holes, 
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where r = (x, y) is the 2D coordinate, m e and rrih are the material-dependent effective masses 
of electrons and holes, assumed isotropic in the x, y plane, E c {v) and E v (r) are the local 
conduction-band and valence-band edges, respectively, and V^(r) is the electrostatic potential 
generated by free carriers and fully-ionized dopants. The above equations are numerically 
integrated to give the set of eigenfunctions (eigenenergies) for electrons and holes e (r) (e e ) 
and 4> h (r) (e h ), respectively. 

After the solution of Eqs. ([1]) and (T5]), we calculate the volumetric total charge density 

p(r) = e(n A (r) - n e (r) + p D (r) - p A (r)), (3) 

where po and pa are the ionized donor and acceptor densities, respectively, and the densities 
of free electrons and holes are 



Here, m e (rrih) is the effective electron (hole) mass along the CSNW axial direction (which 
is in general different from the in-plane mass in Eqs. fTH2l particularly for holes), T is the 
temperature, p is the Fermi level, and = r(fc+i) ^o°° J-^+x * s ^ e complete Fermi-Dirac 

integral of order k, which results from the integration of the parabolic dispersion along the 
growth direction z. 

From the total charge density, we compute the electrostatic potential V solving the Pois- 
son equation with a material-dependent relative dielectric constant s r 

V r [ £r (r)V r y(r)] = -^. (6) 

We use Dirichlet boundary conditions, with the potential on the domain boundaries fixed 
to either an input value representing the voltage of a metallic gate in that position, or zero. 
Several test simulations, with increasing domain size, confirm that the electrostatic potential 
at the CSNW surface is essentially zero when no gate is included. 

The contribution due to exchange and correlation effects is calculated via a local density 
approximation from the densities n e (r) and rih(r) : and added to V. We adopt the well known 



approximate expression reported in Ref. |23| and |2J. In all the simulations presented in the 
following, it turns out to be almost negligible 2 -. For the sake of brevity, here we consider 
the exchange- correlation contribution already included in V. 



The whole procedure is then iterated, with the new potential V inserted in Eqs. ([I]) and 
02]). The self-consistent cycle is stopped when the relative variation of the charge density is 
lower than a given value E p at any position, i.e., 



where A is the area of the simulation domain. 

Once the convergence is achieved, one can calculate the linear charge density (we recall 
that the system is assumed to be translationally invariant along z) by the area integration 



To quantify the charge distribution in the section of the CSNW structure, we also define the 
effective sheet charge density 



where L is one of the two directions shown in Fig.[H which either join two edges or two facets, 
and r(Z) is a point on one of these lines. p s heet represents the sheet charge of an equivalent 
planar heteroj unction, as the charge were uniform along the interface. In our case, p s h ee t is 
in general different in the two directions, and is an indication of the preferential localization 
of the free charge. 

With regards to the numerical approach, the solutions of Schrodinger and Poisson equa- 
tions are obtained through a box integration method^, on a triangular grid with hexagonal 
elements. This choice reproduces the shape of the domain and the symmetry of the (unbi- 
ased) system, and avoids numerical artifacts originated by discretization asymmetries of the 
six domain boundaries, as would be the case, e.g, using a rectangular grid. Formally, the 
partial differential equations Eqs. ([I]), (J2J) and (jfJJ) are integrated on each hexagonal element. 
By applying the divergence theorem, the area integral is converted in a linear integral of the 
flux along the hexagon boundary. A balance between incoming and outgoing fluxes (obtained 
through a first-order finite-differences scheme) of adjacent hexagons connects the unknowns 
on different elements. This results in a symmetric (Hermitian) sparse matrix for the Poisson 
(Schrodinger) equation, whose dimension correspond to the number of hexagons, and with 
six non-zero elements on each row. The matrix for the Schrodinger equation is diagonalized 
through a Lanczos library algorithm 26 while that for the Poisson equation is solved, with 
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the known term obtained from p, via an efficient library routine 27 . The Fermi-Dirac integral 
is computed by using an efficient routine implementing Chebyshev polynomial expansion^. 

Stability of the self-consistent cycle for large and/or heavily doped structures is delicate. 
The usual mixing procedure using an under-relaxation parameter, a, is implemented: at 
iteration i, the electrostatic potential at the current iteration V 1 is substituted by 

V* <- aV i + (l-a)V i -\ (fO) 

In addition, although the Poisson and Schrodinger equations are solved on the same grid 
points, the charge density extends very little into the AlGaAs barriers, due to the large 
depletion operated by surface states. Therefore, the Schrodinger equation need to be solved 
in a much smaller region. Typically we require the wave functions to vanish at the position 
of the doping layer. 

Holes are treated in a single-band approach. As we shall discuss in Sec. IIII C\ the hole gas 
is mostly confined at the GaAs/AlGaAs interfaces, which are directed along the six {110} 
directions. The parabolic dispersion for heterostructures grown along the [110] direction is 
given in22, Eq. 4.4a. The strongly anisotropic mass is heavy in the [110] direction (roughly 
twice the mass along the widely used [100] direction) and light along in-plane direction, here 
corresponding to the vertical growth direction [111], as shown later in Table fl] 

III. RESULTS 

A. Structure and simulation details 

As a prototype structure (see Fig. CP, we simulate a device similar to the one described in 
Ref. y. A GaAs core 80 nm wide is surrounded by a 50 nm wide Alo.3Gao.7As shell, and by 
a 10 nm wide GaAs capping layer. The Alo.3Gao.7As shell is uniformly doped in the center 
with a 10 nm wide layer at a constant doping density pp. Material parameters are shown 
in Tab. [Q p is taken at the mid-gap value of GaAs, and all calculations are performed at 
T = 20K. 

In this Section we report results for two kind of structures, one n-doped and one p-doped. 
In both cases, the occupation probability of minority carriers is negligible, hence, we need to 
solve only the relevant equation, among Eq. ([I]) and ([2]). Schrodinger and Poisson equations 
are solved on a 2D hexagonal domain, discretized with a regular hexagonal tessellation. Our 
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GaAs Alo.3Gao.7As 

E c - E v [eV] 1.43 1.858 

AE C [eV] 0.284 

AE V [eV] 0.144 

m e 0.067 0.092 

m e 0.067 0.092 

m h 0.690 0.731 

m h 0.105 0.124 

e r 13.18 12.24 

TABLE I: Material parameters used in the simulations. AE C and AE V are the conduction and 
valence band offsets, respectively, between the conduction band edge E c and valence band edge E v 
of the GaAs and the Alo.3Gao.7As layers.— m e and m e are the relative electron effective masses 
in the x, y plane, orthogonal to the CSNW axis, and along the axis direction, z, respectively. 
Analogously for and for holes. Note that the electron effective mass is assumed isotropic, 
while hole relative effective masses are strongly anisotropic. 

simulations typically use about 2 x 10 5 hexagonal elements. The Schrodinger equation is 
solved on a sub-set of 5 x 10 4 elements, as explained in the previous section. The under- 
relaxation parameter a ~ 0.05 and the maximum convergence relative error allowed in the 
charge density E p = 10~ 3 . Convergence is typically achieved within 100-^200 iterations 
if the potential V at the first iteration is zeroed. Clearly, the number of iterations can 
be substantially reduced by initializing the potential with a suitable guess, typically, the 
converged potential at a slightly different doping density. Converged results, of course, do 
not depend on this choice. 



B. Electron gas localization 

We now consider a n-doped CSNW, where a free-electron gas is eventually formed in 
the GaAs core or at the inner heterointerface, and specifically investigate the free electron 
gas formation as a function of the doping density pp. Figure [2] (top panel) shows the 
linear free electron density vs po- At a threshold density ~ 1.4 x 10 18 cm -3 the nanowire 
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FIG. 2: Top: Linear charge density of the free-electron gas as a function of the donor density 
Pd- Bottom: sheet charge density of the free-electron gas calculated along the facet-to-facet and 
edge-to-edge directions indicated in Fig. [TJ Arrows indicate three different densities shown in in 
Fig. [3] As the doping density increases, the latter dominates, indicating that the free electrons 
tend to accumulate on the six CSNW edges and form six ID channels. In both graphs, the solid 
(dotted) curves represent the results of self-consistent calculations including (neglecting) exchange- 
correlation contributions. 



starts to be populated, in qualitative agreement with IdJlOl The linear density is a fraction 
of ~ 10 7 cm~ 1 , and increases linearly with doping. Note the tiny difference between the 
calculations with and without the XC potential. In the rest of the paper, all results are 
reported for calculations which include the XC potential. 

A glimpse to the charge density maps reported in Fig. [3] shows that the free electron gas is 
localized in the GaAs core, as expected. Indeed, at sufficiently large doping, the conduction 
band bends down from the surface value until the Fermi level (zero in the right energy scale 
of the right plots of Fig. [3]), which is pinned to the middle of the gap at the CSNW surface, 
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is situated just above the GaAs conduction edge of the core. This is in complete analogy 
to the formation of an inversion layer in a remotely doped planar hetero junction. However, 
here the localization is not homogeneous along the heterointerface, and its pattern strongly 
depends on the doping density, as we discuss below. 

In Fig. [2] (bottom panel) we show the sheet charge density along facets and edges. At 
small doping the two densities are equal, while at larger doping the edge population becomes 
dominant. This evolution is made more clear in the charge density maps of Fig. |3l which 
refer to the three densities labelled (a), (b) and (c) in Fig. |2j At the lowest doping (Fig. Efa), 
when the GaAs core starts to be populated, the charge is distributed deep into the core. 
The distribution is only slightly modulated (right panels) crossing the core along either 
the facet-to-facet of edge-to-edge directions, and slightly depleted in the center of the core. 
Furthermore, the distribution is almost circularly symmetric. Accordingly, the two sheet 
charge densities in Fig. [2] are equal. 

As the doping is increased [Fig. WJd>)] the charge depletion in the center is more pro- 
nounced, and charge moves toward the interfaces. The 2D maps (left panel) shows that the 
distribution is starting to develop a six-fold symmetry, and the sheet charge density along 
the edge-to-edge direction is starting to dominate. For large doping [Fig. Etc)] the charge is 
strongly localized at the edges. This can be seen in the 2D maps (left panel) as well as, more 
quantitatively, in the two profiles (right panels). In this regime the charge density has fully 
developed the six-fold symmetry and resemble a set of coupled quantum wires more than a 
2DEG, since most of the charge is confined in relatively narrow channels at the edges, while 
only a minor part of the charge sits at the facets. 

The charge density is obtained from the occupation of an increasing number of energy 
subbands. The evolution of subband energies with the doping density is shown in Fig.HJwith 
respect to the Fermi energy. As the doping increases, subbands fall more and more below 
the Fermi level. Due to the six-fold symmetry, these states are singly or doubly degenerate, 
as indicated. 

C. Valence band holes 

Nanowires are easily p— doped, for example with Si, which is amphoteric for GaAs.— In 
this section we discuss the formation of the hole gas in the same structure as in the previous 
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FIG. 3: Spatial distribution of the free-electron gas for the three values of the doping density 
indicated as (a), (b), and (c) in Fig. [2] (bottom panel). For each density we show (left) the 2D map 
of the charge, with the GaAs/AlGaAs interface indicated as a solid line, and (right) the charge 
density profile (thin line, left axis, units of 10 15 cm -3 ) and the self-consistent conduction-band 
profile (thick line, right axis, units of eV, \x = 0) along the edge-to-edge direction (right, top) and 
along the facet-to-facet (right, bottom). Note the different scale for different doping density. 

section. We use a parabolic description with the appropriated effective masses, as discussed 
in Sec. [UJ reported in Tab. [U Due to this parabolic model, the difference between electron 
and hole gases to be discussed below reside only in the much heavier effective mass of holes. 

Figure shows the free charge density of holes as a function of the doping level. Similarly 
to the electron gas, there is a threshold doping density 1.4 x 10 18 cm -3 , below which no free 
charge is obtained. Above this threshold the linear charge density (top panel) increases 
linearly, as for conduction electrons. The equivalent sheet charge density (bottom panel) 
along the edges-to-edge and facet-to-facet directions [Fig. [5] of the structure also behave 
similarly to the case of electrons. In particular, at low doping the two densities are similar, 
with the sheet charge density at the edges taking over when the doping increases. 
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FIG. 4: Electron subband energy vs doping density pa- Each bullet represents the subband edge 
resulting from the self-consistent calculation. The lines (thinner for non-degenerate subbands, 
thicker for doubly degenerate subbands) are a guide to the eyes. 

However, it turns out that electrons and holes have a very different behavior and local- 
ization with respect to conduction electrons, particularly at low and intermediate charge 
densities. Figure |6] shows the free hole density profiles at the three densities (a), (b), and (c) 
indicated in Fig. (bottom panel). Although the sheet charge density is nearly equal in the 
two directions at the lowest density, the hole charge density is peaked at the facets, while 
the distribution along the edge-to-edge direction is low at the edges and broader toward 
the core. This is sharp contrast with conduction electrons, where at low doping the charge 
forms a cylindrical distribution in the core. 

At intermediate densities [Fig. E^b)] the gap between the facets is filled, and the hole 
gas is uniformly distributed at the GaAs/AlGaAs interfaces. Only weak modulations are 
present along the interface going around the core in this very interesting intermediate regime. 
Therefore, in this case the hole gas indeed resembles a 2D hole gas which is six-fold bent 
around the wire. This regime has no counterpart for conduction electrons. At the largest 
density [Fig. Etc)], instead, the charge is concentrated at the edges, similarly to the case of 
electrons. 

The evolution of the hole subbands with doping density is shown in Fig. [7J The behavior 
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FIG. 5: Top: Linear charge density of the free-hole gas as a function of the acceptor density pa- 
Bottom: sheet density of the free-hole gas calculated along the facet-to-facet and edge-to-edge 
directions indicated in Fig. [TJ Arrows indicate three different densities shown in in Fig. 

is clearly similar to that of conduction electrons, except for the large density of levels due to 
the larger mass. Two sets of levels can be clearly recognized, corresponding to levels with 
one additional nodal surface in the radial direction. 

D. Origin of localization patterns 

To rationalize the behavior of electrons and holes vs doping, particularly at low density, 
we can consider the following energy contributions: 1) repulsive Coulomb energy between like 
particles, which requires maximization of the average inter-particle distance and therefore 
favors localization at the edges of the inner GaAs/AlGaAs interface, 2) vertical localization 
energy, that is the kinetic energy required to localize charge at the GaAs/AlGaAs interfaces, 
analogously to a planar heterojunction, and 3) lateral localization energy, that is the energy 
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position (nm) 

FIG. 6: Spatial distribution of the free-hole gas for the three values of the doping density indicated 
as (a), (b), and (c) in Fig. [5] (bottom panel). For each density we show (left) the 2D map of the 
charge, with the GaAs/AlGaAs interface indicated as a solid line, and (right) the charge density 
profile (thin line, left axis, units of 10 15 cm -3 ) and the self-consistent valence-band profile (thick 
line, right axis, units of eV, /j, = 0) along the edge-to-edge direction (right, top) and along the 
facet-to-facet (right, bottom). Note the at the lowest doping (a) the free holes are localized mainly 
along the facets, contrary to the electron case, and at intermediate doping (b) the free holes form 
an hexagonal ring, a regime which have no counterpart for conduction electrons. 

required to localize parallel to the interface, along a facet or at an edge, the latter being 
clearly larger. 

At low density, when the Coulomb contribution is small, conduction electrons distribute 
over the core, since vertical localization energy is large. Lateral localization would be also 
dominant over the small Coulomb energy in this regime, hence the isotropic, cylindrical sym- 
metry. The charge gets localized near the GaAs/AlGaAs interfaces only if the density and, 
hence, the average Coulomb energy is so large that localization at the edges is energetically 
favorable. In this regime, indeed, the small additional cost in lateral localization energy 
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FIG. 7: Hole subband energy against doping density pa- Each bullet represents the subband 
threshold resulting from a self-consistent calculation. The lines are a guide to the eyes. 

at the edges with respect to facets is overcompensated by gain in Coulomb energy due to 
inter-edge repulsion. For holes, instead, the small kinetic energy (one order of magnitude 
smaller than for conduction electrons) makes it easier to localize at the accumulation layer 
near the interface even at very low densities. Since the Coulomb energy is anyway small 
in this regime, at the interfaces hole density may localize at the facets, rather than at the 
edges, where the lateral confinement energy would be larger. 

We finally note that the localization depends essentially on the amount of free charge, 
which in the present calculation is tuned by the doping density. Alternatively, one could fix 
the doping density to a large value, and deplete the electron (hole) gas using an external 
gate grown around the wire, recovering the different localization regimes shown above using 
a negative (positive) voltage applied to it. We have explicitly calculated the charge evolution 
by simulating an all around gate, and indeed recovered very similar results (not shown). In 
the next chapter, we will show instead what happens by application of a back-gate to a NW. 



15 



E. Biased system 

Transport properties of NWs can be measured by integrating them in planar field effect 
transistors (FET)^^. In such devices a NW is placed on a substrate formed from a dielectric 
layer which separates the wire from another layer acting as a global back gate. The applied 
gate voltage depletes or populates the NWs depending on the sign, and hence, it can be used 
to modulate the conductance between the source and drain electrodes. We next simulate the 
effect of such a back gate assuming that it is in remote contact with one of the external facets 
of the NW, hereafter referred to as the bottom facet (see inset in Fig. [9]). The potential 
is set to V = —eV g on the basis edge, where V g is the applied electrode voltage in the 
range V g — ±1 V. Since the boundary of our simulation domain corresponds to the CSNW 
surface and it is at least 60 nm distant from the free carriers, we assume that the possible 
modulations of V along the edges have little influence on the carrier localization in the core. 

In Fig. [8] we show the electron charge density profiles at four selected gate voltages for 
two NWs with a low (left) and high (right) density of dopants. As anticipated, the charge 
density no longer shows the symmetric distributions obtained at zero voltage, and tends to 
accumulate at the top and bottom of the wire section for negative and positive voltages 
respectively. Interestingly, for the NW with higher density of dopants we observe how the 
six quasi-lD channels obtained at zero voltage can be switched to four and two channels by 
tuning the voltage to V g = —0.1 V and V g = ±0.4 V. Also note the different localization for 
negative voltages at low or high doping, where charge is distributed in a broad channel or 
in two narrower channels, respectively. 

The quantitative effect of the gate voltage can be observed in Fig. where we represent 
the linear charge density as obtained by integrating the volume charge density over the 
three domains depicted in the inset, as a function of the applied voltage. The total linear 
charge density is also plotted (black solid line). Clearly, the three curves coincide at zero 
bias due to symmetry. A positive gate voltage does not deplete the top regions of the wire 
in favor of a larger electronic concentration in the bottom of the wire. Instead, the charge 
density raises quasi-linearly in all regions, although this increase is much pronounced in the 
bottom regions Rl and R2, while in region R3 the charge remains nearly constant at the 
zero bias value. On the other hand, the evolution of the charge density is not linear when 
the NW becomes depleted by applying negative voltages. This fact can be associated with 
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FIG. 8: Spatial distribution of the free-electron gas in two CSNWs with doping densities 1.4 x 
10 18 cm -3 (left) and 1.64 x 10 18 cm -3 (right) at four selected values of the gate voltage. Depending 
on the sign of the voltage, the free electrons tend to localize near the top/bottom parts of the wire. 

the smaller density of subbands to populate at such low-density regimes. Finally, we notice 
that, despite the different density of dopants, at high positive voltages both systems reach 
the same density regime with the electrons concentrated at the two bottom edges of the wire 
(see, e.g., bottom plots in Fig [8]). 

Figures [10] and [11] show the corresponding results for hole populated NWs with two 
different densities of dopants. Clearly, the applied voltage induces the opposite effect than 
in the electronic system, increasing (decreasing) and localizing the charge density at the 
bottom (top) of the NW at negative (positive) voltages. In parallel with the result obtained 
at zero voltage, the NW with higher density of p-dopants shows qualitatively the same 
behavior as its n-doped counterpart. Thus, in the right panels of Fig. [10] it can be seen 
that, again, a NW with two and four quasi- ID hole channels can be attained by tuning the 
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FIG. 9: Linear electron density integrated over the three domains of the NW illustrated in the 
inset, as a function of the gate voltage. Solid blue lines are used for the domain labeled as Rl, 
dashed red lines for R2, and dotted yellow lines for R3. The total electron density is also plotted 
in thick solid black lines. Two densities of n-dopants are represented (a) 1.4 x 10 18 cm -3 and (b) 
1.64 x 10 18 cm- 3 . 



applied voltage. Nevertheless, the trend of the hole density to localize in the facets at low 
density regimes enables the observation of new voltage-induced charge arrangements when 
the density of dopants is lower. For instance, at V g = 0.05 V one can see that the charge 
density is concentrated in the center of the three top facets of the NW, while at V g = —0.3 V 
a 2D hole gas is clearly formed over the bottom facet. 

To conclude, in Fig. [11] we show the effect of the gate voltage on the linear hole density 
integrated over the three regions of the NW. The results are equivalent to the electronic 
case, with the role of the three regions reversed due to the opposite sign of the carriers, and 
can be rationalized along the same lines. 
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FIG. 10: Spatial distribution of the free-hole gas in two CSNWs with doping densities 1.45 x 
10 18 cm -3 (left) and 1.65 x 10 18 cm -3 (right) at four selected values of the gate voltage. New charge 
arrangements are observed not present in the electron-gas case, e.g. the charge accumulation at 
the three top facets of the wire (V g = 0.05 V). 

IV. CONCLUSIONS 

We have determined the properties of an electron/hole gas in a realistic model of radial 
hetero junctions formed in a modulation doped GaAs/AlGaAs CSNWs, emphasizing the 
different localization and symmetry which can be realized as a function of the doping level 
and/or by a back-gate. Contrary to planar heteroj unctions, conduction electrons do not form 
a 2DEG well localized at the GaAs/AlGaAs interface, but rather show transitions between 
different charge distributions across the section of the CSNWs. For conduction electrons, 
these are similar to phases predicted in GaN/AlGaN CSNWs 2 =. For p-doped structures, the 
heavier effective mass result in additional phases for the hole gas which have no counterpart 
for conduction electrons. In particular, it is possible to form a nearly uniform six-fold bent 
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FIG. 11: Linear hole density integrated over the three domains of the NW illustrated in the inset, 
as a function of the gate voltage. Solid blue lines are used for the domain labeled as Rl, dashed red 
lines for R2, and dotted yellow lines for R3. The total hole density is also plotted in thick solid black 
lines. Two densities of p-dopants are represented (a) 1.45 x 10 18 cm" 3 and (b) 1.65 x 10 18 cm -3 . 

2D hole gas at the GaAs/AlGaAs interfaces. 

In general, however, our results show that the electron and hole gases in these structures 
deviate substantially from any isotropic distribution at the GaAs/AlGaAs heterointerfaces. 
It is particularly important that, depending on the density and possibly the applied gate 
potential, the charge density reshape between quasi-lD and quasi-2D distributions. This 
is likely to be reflected in multiple and tunable time-scales of scattering mechanisms.— In 
most cases, therefore, we expect complex behavior in possible future transport experiments. 
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